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Ground state properties of proton-rich oAd-Z nuclei in the region 55 < Z < 73 are 
studied in the relativistic mean field (RMF) theory. The RMF equations are solved by using 
QQ ■ the expansion method in the Harmonic-Oscillator basis. In the particle-particle channel, we 

use the state-dependent BCS method with a zero-range 5-force, which has been proved to 
be effective even for neutron-rich nuclei. All the ground state properties, including the one- 
. proton separation energies, the ground state deformations, the last occupied proton orbits 

^ ' and the locations of proton drip line, are calculated. Good agreement with both the available 

\^ ' experimental data and the predictions of the RHB method are obtained. 
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The structure and decay modes of nuclei at and beyond the proton drip hne 
represent one of the most active areas in both experimental and theoretical stud- 
ies of exotic nuclei with extreme isospin ratios. Recent reviews can be found in 
Refs.^-*'^^ and references therein. For each element, the number of neutrons that can 
be supported by the nucleus is limited if the nucleus is to remain stable with respect 
Pi , to nucleonic emission. These limits define the proton and neutron drip lines. The 

main difference between the proton and the neutron is the existence of the Coulomb 
potential for the former. Consequently, the proton drip line lies much closer to the 
stability line than the neutron drip line. For example, the lightest particle-stable 
\ neon isotope is ];QNe7^'' , while as heavy an isotope as fo^e24^^ is known to exist. 

The existence of the Coulomb barrier also traps the wave function of the proton 
in the nuclear region, which explains why we have not yet clearly observed halo 
phenomena in proton-rich nuclei. 

The proton drip line has been almost fully mapped up to Z = 21^^ experimen- 
tally. No proton-unbound nuclei have been directly observed in this region, which is 
not surprising, since the Coulomb barrier is relatively low. It has been argued that 
ground state proton decay can only be detected directly for nuclei with Z > 50, where 
the relatively high potential energy barrier causes nuclei to survive long enough to be 
detected experimentally. Detailed studies of ground state proton radioactivity have 
been reported for odd-Z nuclei mainly in the two mass regions, 51 < Z < 55 and 



*•* E-mail: lsgengO@rcnp.osaka-u.ac.jp 
**-* E-mail: toki@rcnp.osaka-u.ac.jp 
**•) E-mail: mengj@pku.edu.cn 



typeset using PT^TeX.cIs (Ver.0.89) 



2 



Q9 < Z < 83. The experimental features observed for most of these nuclei have been 
explained by assuming them to have spherical shapes.^) The half-lives for proton de- 
cay were then evaluated using the semi-classical WKB method or standard reaction 
theory within the distorted wave approximation. However, to obtain proton decay 
rates in the transitional region beyond the Z = 50 shell closure (^"^^I and ^^^'^^"^Cs) 
and in the region of light rare- earth nuclei (^^^Eu and ^^O'^^^Ho) requires calcula- 
tions assuming significant prolate deformations. Recently, an exact formalism was 
developed in Refs.,'^^'^^ where the half-lives are evaluated with the assumption that 
the emitted proton exists in a deformed single particle Nilsson level. This formalism 
was applied^)"^) to analyze all data presently available for odd proton emitters with 
an even number of neutrons and the known odd-odd proton emitters ^-^^Cs, -"^^^Ho 
and ^"'''^Lu.'^) This formalism, however, does not predict proton separation energies, 
i.e. the model does not predict which nuclei arc likely to be proton emitters. 

In the past decade, the relativistic mean field (RMF) theory has been successfully 
applied to the study of nuclear properties throughout the Periodic Table. General 
reviews of the RMF theory and its applications in nuclear physics can be found in 
Refs.^)"^^) To study exotic nuclei whose Fermi surfaces are close to the threshold, 
several issues, including the continuum and the pairing correlation, must be treated 
carefully. ^^^"^'^^ Luckily, in the case of proton-rich nuclei, because the nuclear wave 
function is well confined in the nuclear region due to the existence of the Coulomb 
potential, one only needs to take care of the pairing correlation. The use of a zero- 
range (5-force in the particle-particle channel, V = — Vo<5(ri — f-z), has proved to be an 
efficient and economical way to accomplish this.^^^"^*^) (For detailed discussions of 
more general cases, including neutron-rich nuclei, we refer the reader to Refs. 1^)^20) 
and references therein.) In our present work, the deformed RMF-I-BCS method 
with a zero-range 5-force in the pairing channel^'^^"^^^ is applied to the analysis of 
deformed odd-Z proton-rich nuclei with 55 < Z < 73. In the mean-field part, two of 
the most successful parameter sets, TMA^^) and NL3,^^) are used. 

The RMF theory was first applied to study deformed proton emitters by Lalazis- 
sis et al. In Refs.,^^)"^^) the proton drip-line nuclei with 31 < Z < 49 and 51 < Z < 
73 are studied within the framework of the Relativistic Hartree-Bogoliubov (RHB) 
model. In the mean-field part, the NL3^"^^ parameter set is used. In the particle- 
particle channel, the pairing part of the Gogny force with the DIS set^^^ has been 
adopted. The calculated one-proton separation energy and other relevant proper- 
ties are in good agreement with the available experimental data. Such a success 
shows that the RMF theory could be an appropriate method even for the descrip- 
tion of nuclei with extreme isospin ratios. This observation also motivates us to 
study deformed proton emitters within a different framework, that of the deformed 
RMF-I-BCS method, due to the following considerations. First, we would like to see 
whether such a success is dependent on the choice of a particular combination of 
parameter sets, i.e. NL3-I-D1S. Second, we want to see whether the argument that 
a zero-range (5-force is an efficient and economical interaction to treat the pairing 
correlation holds in such extreme cases, i.e. at and/or near the proton drip line. 
Finally, many new theoretical and experimental works, including those reported in 
Refs.,^^'^)'^'^)'^^) to name just a few, have been carried out since the last series of 



3 

works^^^"^^-* in the relativistic mean field theory. Because in general, the BCS method 
is thought of as an approximation to the full Bogoliubov transformation, we would 
like to compare our results with the existing RHB results. 

The article is organized as follows. In §2, we give a brief introduction to the 
deformed RMF+BCS method. In §3, we present our results for odd-Z proton-rich 
nuclei with 55 < Z < 73 and compare our results with the available experimental 
data and those of the RHB method. §4 is devoted to the summary of this paper. 

§2. The deformed RMF+BCS method 

Our RMF calculations have been carried out using the model Lagrangian density 
with nonlinear terms both for the a and lu mesons, as described in detail in Ref.^^^ 
The Lagrangian density is given by 



—ni^^n/^" + ^m^w^o;^ + ^g^iujf.ij'^f - 5a;V'7''V'^M 

-\f,,F^^ - eViTM^^^'^V', (2-1) 

where the field tensors of the vector mesons and of the electromagnetic field take 
the forms 

= d^ujy — dijUjfj,, 
^M-= d^pi-d,p';^-2g,e-'^-py,, (2-2) 
Fjxv — d^Ai, — diiA^, 

and other symbols have their usual meanings. Based on the single-particle spec- 
tra calculated with the RMF method, we perform a state- dependent BCS calcula- 
tion.^^)' The gap equation has a standard form for all the single particle states, 

= _1 V (2-3) 

where e^i is the single particle energy and A is the Fermi energy, whereas the particle 
number condition is given by 2 ^ = AT. In the present work, for the pairing 

A;>0 

interaction we use a zero-range 5-force, 

V = -Vo5{ri-r2), (2-4) 

with the same strength Vq for both protons and neutrons. The pairing matrix element 
for the zero-range (5-force is given by 

Vij = {fi\V\jj) - {fi\V\jj) = -Vo J d^r - ^M'i'-M ' (2-5) 
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with the nucleon wave function in the form 

A detailed description of the deformed RMF+BCS method can be found in Ref.^'^) 
In the present study, nuclei with both even and odd numbers of protons (neu- 
trons) need to be calculated. We adopt a simple blocking method, in which the 
ground state of an odd system is described by the wave function 

a^JBCS) = al J] (^'i + Vk44')\''^)- (2-7) 

Here, |vac) denotes the vacuum state. The unpaired particle sits in the level ki 
and blocks this level from pairing correlations. The Pauli principle prevents this 
level from participating in the scattering process of nucleons caused by the pairing 
correlations. As described in Ref.,^^) in the calculation of the gap, one level is 
"blocked", 

A, = -'- ^ , ■ (2-8) 

fcVfei>o y (£fe' - A)2 + Al, 

The level ki has to be excluded from the sum, because it cannot contribute to the 
pairing energy. The corresponding chemical potential. A, is determined by 

Ar = l + 2 ^fc- (2-9) 

This blocking procedure is performed at each step of the self-consistent iteration. 

§3. Deformed ground-state proton emitters with 55 < Z < 73 

In the present work, we use the mass-dependent parameter set TMA^^) for the 
RMF Lagrangian. Results with the parameter set NL3^'*^ are also presented for 
comparison. For each nucleus, first the quadrupolc constrained calculation'^"^) '^^^ is 
performed to obtain all possible ground state configurations, and then we perform 
the non-constrained calculation using the quadrupole deformation parameter of the 
deepest minimum of the energy curve of each nucleus as the deformation parameter 
for our Harmonic-Oscillator basis. In the case that there arc several similar minima 
obtained from the constrained calculation, we repeat the above-described procedure 
to obtain the configuration with the lowest energy as our final result. The calculations 
for the present analysis have been performed using an expansion in 14 oscillator shells 
for fermion fields and 20 shells for boson fields for isotopes with Z > 70, i.e. Ta and 
Lu isotopes. While for isotopes with Z < 70, 12 shells for fermion fields and 20 
shells for boson fields are used. In the latter case, we use 12 shells for fermions to 
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save computation time. Convergence has been tested for all the quantities calculated 
here. Following Ref.,^^^ we fix Hujq = 41 A~^/^ for fermions. 

Whenever the zero-range 6 force is used in the BCS framework, a cutoff proce- 
dure must be applied, i.e. the space of the single-particle states in which the pairing 
interaction acts must be truncated. This is not only to simplify the numerical cal- 
culation but also to simulate the finite-range (more precisely, long-range) nature of 
the pairing interaction in a phenomenological way.^^^'^^^ In the present work, the 
single-particle states subject to the pairing interaction are confined to the region 
satisfying 

ei-X< E,^t, (3-1) 

where Cj is the single-particle energy, A the Fermi energy, and E'cut = 8.0 MeV. We 
find that increasing Ecut from 8.0 MeV up to 16.0 MeV. followed by a readjustment 
of the pairing strength Vq, docs not change the results, and therefore none of our 
conclusions will change. The same pairing strength is used for both protons and 
neutrons and is fixed by requiring that the experimental one-proton separation en- 
ergies of several selected nuclei be reproduced by our calculations, including "^^'''Ta, 
-'^^''Tm, and ^-"^^Cs. A slight readjustment of the pairing strength, say %5 , does not 
change our conclusion here. More specifically, for calculations with the NL3 set, we 
take Vq = 393.0 MeV fm^ for Z>71, and Vq = 420.0 MeV fm^ otherwise. For calcu- 
lations with the TMA set, we take Vq = 393.0 MeV fm^ for Z > 71, and Vq = 375.0 
MeV fm^ otherwise. We should mention that the small difference between Vq for 
isotopes with Z > 71 and isotopes with Z < 71 is partly due to the different number 
of shells used in the expansion method. 

In the process of proton emission, the valence proton tunnels through the Coulomb 
and centrifugal barriers, and the decay probability depends strongly on the energy of 
the unbound proton and on its angular momentum. In rare-earth nuclei, the decay 
of the ground state by direct proton emission competes with (3'^ decay; for heavy nu- 
clei, fission or a-dccay can also be favored. In general, ground state proton emission 
is not observed immediately beyond the proton drip line. For small values of the 
proton separation energy, the width is dominated by /3+ decay. On the other hand, 
large separation energies result in extremely short proton emission half-lives, which 
are difficult to observe experimentally. For a typical rare-earth nucleus, the separa- 
tion energy window, in which ground state proton decay can be directly observed, 
is about 0.8-1.7 MeV.^) 

The Nilsson quantum numbers of the last occupied proton orbit are taken to be 
the same as the dominant component in the expansion of this wave function in terms 
of the axially-symmetric Harmonic-Oscillator basis. ^^-^ In theoretical calculations, 
the spectroscopic factor of the corresponding odd-proton orbit, k, is defined as the 
probability that this orbit is found empty in the daughter nucleus with one less 
proton. If the difference between the ground state configurations of the mother 
and the daughter nuclei is ignored, the spectroscopic factor for level k is calculated 
as 

Sk = maMl'' = 4- (3-2) 
Here, is related to the occupation probability v'^ for level k in the daughter nucleus 
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through the weh-known relation in the BCS formalism 

ul + vl = l. (3-3) 

In what follows, we discuss the details of our numerical results. 

3.1. Lutetium (Z=71) and Tantalum (Z=73) 

The one-proton separation energies for Lu and Ta isotopes are displayed in Fig. 

1 as functions of the mass number A. The drip-linc nucleus for the Ta isotopic chain 
is predcitcd to be -'^^'^Ta by NL3 and ^^-"^Ta by TMA. While for the Lu isotopic chain, 
both TMA and NL3 predict the drip-line nucleus to be ^^^Lu. We also compare the 
calculated separation energies with the available experimental transition energies 
for ground state proton emission in ^^"Lu, ^^^Lu,^^") ^^^Ta,*^^) ^^^Ta'^'^'> and ^^^Ta.^^^^ 
The corresponding RHB results are taken from Ref.^^^ In Table I, we list the ground 
state properties of Lu and Ta isotopes. The separation energy window, 0.8-1.7 MeV, 
extends to include those nuclei for which a direct observation of ground state proton 
emission is in principle possible on the basis of our calculated separation energies. 
Good agreement between our calculations and both the experimental data and the 
RHB results is clearly seen. However, for the odd-odd nuclei ^^°Lu and ^^^Ta, our 
calculations predict values of Sp smaller than the experimental values. This may be 
explained by the possible existence of a residual proton-neutron pairing interaction. 
In both nuclei, an extra positive energy of about 0.3 MeV is needed to increase the 
calculated Sp value to the experimental value. We find that if we reduce the pairing 
strength for odd-odd nuclei by about 5-10% compared with odd-even or even-even 
nuclei, wc can fit the experimental transition energy quite well. Wc sec below that 
the same conclusion also holds for other odd-odd nuclei. Although it has long been 
known that the proton-neutron pairing is important for the proper description of 
nuclei,^^)'^^) it was believed that it only affects those proton-rich nuclei with N ^ Z . 
However, recently Simkovic et al.^^^ showed that even for nuclei with N — Z = 12, 
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Fig. 1. The one-proton separation energies for Lu {Z=71) and Ta {Z=73) isotopes at and beyond 
the proton drip line. The experimental value for the proton separation energy corresponds to 
the opposite of the ground state transition energy, Ep. 
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Table I. Lu {Z=71) and Ta (Z=73) ground state proton emitters. The results of RMF+BCS/TMA 
calculations (third to sixth columns) for the one-proton separation energy, Sp, the mass 
quadrupolc deformation, P2, and the deformed single-particle orbit occupied by the odd va- 
lence proton, are compared with the predictions of the RHB model^^^'^^^ (seventh to tenth 
columns) and with the experimental values, Sp = —Ep, where Ep is the experimental transition 
energy. The RMF-I-BCS/TMA and RHB spectroscopic factors are displayed in the sixth and 
tenth columns, respectively. All energies are in units of MeV. 





RMF-hBCS/TMA 


RHB 


Exp. 


iV 


Sp 


P2 


p orbital 




Sp 


P2 


p orbital 




Sp 


^*^Lu 


78 


-1.946 


-0.166 


5/2^ [512] 


0.574 


-1.77 


-0.158 


7/2" [523] 


0.60 




isoLu 


79 


-1.613 


-0.129 


5/2" [523] 


0.549 


-1.31 


-0.153 


7/2- [523] 


0.61 


-1.261(4)''°) 


is^Lu 


80 


-1.239 


-0.119 


5/2" [532] 


0.541 


-1.24 


-0.151 


7/2" [523] 


0.58 


-1.233(3)'*°) 


i^^Lu 


81 


-1.141 


-0.053 


5/2" [532] 


0.496 












i^^Lu 


82 


-1.156 


0.002 


7/2" [523] 


0.463 












i^^Lu 


83 


-0.720 


-0.050 


5/2" [532] 


0.513 












is^Ta 


82 


-1.698 


0.007 


9/2" [514] 


0.374 


-1.677 


0.000 


11/2" 


0.60 


-1.765(10)''i' 


is^Ta 


83 


-1.462 


-0.048 


3/2" [521] 


0.437 


-1.129 


0.000 


3/2+ 


0.51 


-1.007(5)*"' 


i^^Ta 


84 


-0.974 


0.036 


9/2" [514] 


0.381 


-1.040 


0.000 


11/2" 


0.42 


-0.927(7)*^) 


i^^Ta 


85 


-0.725 


0.077 


9/2" [514] 


0.470 













the proton-neutron pairing cannot be ignored. 

With regard to the deformation, both our RMF-j-BCS method and the RHB 
model predict oblate shapes for Lu proton emitters and similar values for the ground 
state quadrupole deformations. Recent calculations by Fcrrcira ct al.^^ show that the 
proton decay in ^^^Lu can be described very well as a decay from a, K = b/2~ ground 
state with an oblate deformation for which —0.18 < /32 < —0.14. In another work,^) 
the author finds that the proton decay in "^^'^Lu is described very well as a decay from 
a, K = 5/2^ ground state with an oblate deformation satisfying —0.17 < (32 < —0.16. 
These calculations are in reasonable agreement with our predictions (see Table I). 
While the RHB method^^^ predicts that the last occupied orbit is 7/2~[523] for both 
nuclei. If we look into the single-particle spectra of these two nuclei more closely, 
we find that in addition to the 5/2~ possibilities, which are listed in Table I, the 
orbits 1/2+ [400] and 7/2" [523] are also possibly blocked. The total energy difference 
between these blocking choices are £^5/2 — -£^1/2 = 0.061 MeV and £'5/2 — -E7/2 = 0.303 
MeV for ^^^Lu; £5/2 - £1/2 = -0.040 MeV and £5/2 - £7/2 = 0.242 MeV for ^^^Lu, 
where E is the binding energy for our calculations with the TMA parameter set. 

In calculations using the RHB model, ^^-^ Ta isotopes are assumed to be spherical, 
while they are slightly deformed in our calculations. At first sight, the RMF-I-BCS/TMA 
and RHB models seem to predict quite different p-orbitals for ^^^Ta, ^^^Ta, and 
^^^Ta, but in fact these states are quite close in our calculations. The theoreti- 
cal predictions are compared with the corresponding experimental assignments for 
the ground state configuration and spectroscopic factor. The experimental data for 
^^^Ta,'^^) i^^Ta^) and i^'^Ta^) are, respectively, the p-orbital 11/2" with 0.58(20), 
the p-orbital 3/2+ with 0.67(16) and the p-orbital 1/2+ with 0.74(34). We notice 
that the experimental assignments are somewhat different both from our predictions, 
in which 155-I57rpjj 

are slightly deformed, and from the RHB calculations, in which 
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Fig. 2. The same as Fig. d but for Ho {Z=67) and Tm (Z=69) isotopes. 



spherical shapes are assumed. There are two possible explanations for this disagree- 
ment. First, these states are close to each other in our calculations and therefore 
it is difficult to make a clear distinction between them. Second, deformation could 
slightly change the state occupied by the odd-proton. Because the nuclei ^^^~^^^Ta 
are not strongly deformed, a different combination of deformation and the occu- 
pied state can give the same experimentally observed values, i.e. the experimental 
transition energies and decay half- lives. 

3.2. Holmium (Z=67) and Thulium (Z=69) 

In Fig. 2, we plot the one-proton separation energies for Ho and Tm isotopes 
as functions of the mass number A. The available experimental data for ^'^^Tm,^''^ 
i46rpjj^^47) i45rp^^48) mojjq^s) and ^^^Ho,^^^ together with the results of the RHB 

calculations, are also shown for comparison. The predicted drip-line nucleus for 
the Tm isotopic chain is ^^^Tm, while for the Ho isotopic chain, NL3 predicts ^^^Ho 
and TMA predicts ^^''Ho. From Fig. [2 we can see that both RMF-hBCS/TMA and 
RHB predict a smaller transition energy for ^^^Ho, while RMF-I-BCS/NL3 predicts 
a value of Sp that is closer to the experimental value. Nevertheless, all theoretical 
calculations fail to predict the increase of the one-proton separation energy from 
^^^Ho to ^^'^Ho. This can be corrected by reducing the pairing strength by a few 
percent, as discussed above. However, this is not our aim here. It is quite clear 



Table II. The same as Table □ but for Ho (Z=67) and Tm (Z=69) isotopes. 





RMF+BCS/TMA 


RHB 


Exp. 


iV 


Sp (32 p orbital 


Sp P2 p orbital u'^ 


Sp 


"°Ho 73 
"^Ho 74 


-1.223 0.341 7/2" [523] 0.592 
-0.885 0.339 7/2" [523] 0.594 


-1.10 0.31 7/2" [523] 0.61 
-0.90 0.32 7/2" [523] 0.64 


-1.086(10)"^' 
-1.169(8)-™> 


"^Tm 76 
"f^Tm 77 
"^Tm 78 


-1.458 -0.220 7/2" [523] 0.458 
-1.238 -0.207 7/2" [523] 0.460 
-0.958 -0.186 7/2" [523] 0.467 


-1.43 0.23 7/2" [523] 0.47 
-1.20 -0.21 7/2" [523] 0.50 
-0.96 -0.19 7/2" [523] 0.55 


-1.728(10)''®' 
-1.120(10)'"') 
-1.054(19)''°' 
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that, even after changing the pairing strength, we cannot account for the trend from 
^^^Ho to ^^'^Ho, which has been attributed to the lack of a residual proton-neutron 
interaction in our above analysis. In Table II, we list the properties of the Ho and 
Tm proton emitters. We note that our RMF+BCS/TMA calculations predict an 
oblate shape for ^^^Tm, while the RHB method predicts a prolate shape. In fact, in 
our constrained calculation, both prolate and oblate shapes are possible for ^^^Tm 
and ^^^Er. This is illustrated in Fig. 3, where we plot the binding energy per 
particle for ^^^Tm and -""^^Er as functions of the quadrupolc deformation parameter, 
/32m.- The binding energies are obtained from the RMF+BCS/TMA calculations 
performed by imposing a quadratic constraint on the quadrupole moment. We find 
two similar minima around (^2 ~ —0.22 and P2 ~ 0.25 for both ^^^Tm and ^^^Er. 
More specifically, for ^^^Tm the two minima are £^7^4 = —7.906 at /?2m = 0.287 
and E/A = -7.903 at P2m = -0.211, while for ^^''Er they are E/A = -7.767 at 
P2m = 0.303 and E/A = —7.768 at /32m = —0.196. The reason that we choose the 
oblate shape instead of the prolate shape for ^^^Tm is two-fold. First, we would like 
to choose the same shape for both ^^^Tm and ^^"^Er. Second, the oblate shape for 
is more consistent with its neighbors (also see Fig. 7). We note that, unlike 
Lu and Ta isotopes, here both RMF+BCS/TMA and RHB predict the same last 
occupied odd-proton state for Ho and Tm isotopes (also see Table II). 
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Fig. 4. The same as Fig. □ but for Eu (Z=63) and Tb (Z=65) isotopes. 



3.3. Europium (Z=63) and Terbium (Z=65) 

In Fig. 4, we plot the one-proton separation energies for Tb and Eu isotopes. The 
results of the RHB modeP^^ and the only available experimental data for ^^ig^so) 
are also shown. The drip-line nucleus for the Tb isotopic chain is predicted to be 
^"^"Tb by NL3 and to be ^^^Tb by TMA. While for the Eu isotopic chain, the drip- 
line nucleus is predicted to be ^^^Eu by TMA and ^^^Eu by NL3. In this region, only 
one proton emitter, ^^^Eu, has been reported. However, based on our calculations, 
there are three other possible proton emitters, ^^''Eu, ^^^Tb and ^^^Tb. We list 
their properties in Table ITTII We notice that the experimental transition energy for 
i^^Eu is 0.950(8) MeV, while our RMF-hBCS/TMA calculations predict 0.498 MeV. 
This can be corrected by adjusting the pairing strength slightly. Considering this 
correction, we see that RHB and RMF-I-BCS predict similar proton emitters for 
Tb and Eu isotopes. As we have seen from the results presented to this point, the 
description of proton emitters is quite sensitive to the pairing strength. We hope 
that further experimental measurements of proton emitters in the entire region can 
help us understand more about the pairing interaction in exotic nuclei. 



Table III. The same as Table |I| but for Eu (^=63) and Tb (Z=65) isotopes. 





RMF+BCS/TMA 


RHB 


Exp. 


N 


Sp (32 p orbital u'^ 


Sp (32 p orbital u'^ 


Sp 


^™Eu 67 
Eu 78 


-0.778 0.420 5/2+ [413] 0.464 
-0.498 0.411 5/2+ [413] 0.323 


-1.22 0.34 5/2" [532] 0.44 
-0.90 0.35 5/2+ [413] 0.44 


-0.950(8)'^°' 


135 rp^^ 70 

136 rpj^ 71 


-1.010 0.365 3/2+ [411] 0.879 
-0.653 0.374 3/2+ [411] 0.862 


-1.15 0.34 3/2+[411] 0.62 
-0.90 0.32 3/2+ [411] 0.65 





3.4. Praseodymium (Z=59) and Promethium (Z=61) 

In Fig. 5, the one-proton separation energies for Pm and Pr isotopes are plotted 
as functions of the mass number A, together with the results of the RHB model. 
The drip-line nuclei are predicted to be ^^^Pm and ^^Spj, i^q^j^ TIMA and NL3. 
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Fig. 5. The same as Fig. Q but for Pr {Z—59) and Pm {Z=61) isotopes. 
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Fig. 6. The same as Fig. but for Cs (Z=55) and La {Z—57) isotopes. 



No experimental proton emitters have been reported for these two isotopic chains. 
Based on our calculations, possible proton emitters in this region are ns-i^Opj. ^^^^ 
I24,i25pj^_ The properties of these nuclei are listed in Table IV. 

3.5. Cesium (Z=55) and Lanthanum (Z=57) 

In Fig. El we plot the one-proton separation energies for Cs and La isotopes. 
The results of the RHB model^'^) are also shown for comparison. The drip-line nuclei 
are predicted to be ^^^Cs and ^^^La by TMA, ^^^Cs and ^^^La by NL3. In Fig. 6, 
we see that the latest reported proton emitter, ^^^La,^°^'^^^ is reproduced quite well 
by RMF-FBCS/NL3 calculations. Both RHB and RMF-^BCS fail to reproduce the 
experimental data for ^^^Cs. As mentioned above, ^^^Cs has an odd number of pro- 
tons and an odd number of neutrons. Since — Z is only 2, we expect a relatively 
strong interaction between the odd proton and the odd neutron. Compared with an 
odd-even system, the additional interaction increases the binding energy somewhat. 
Because the two mean-field models, RHB and RMF-I-BCS, do not include any resid- 
ual proton- neutron interaction, they cannot reproduce the inversion of separation 
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Table IV. The same as Table HI but for Pr {Z=59) and Pm (Z=61) isotopes. 





RMF+BCS/TMA 


RHB 


Exp. 


N 


Sp /32 p orbital u'' 


Sp (32 p orbital u'' 


Sp 


^^'■*Pr 60 
i20pr 61 


"1.354 0.381 3/2" [541] 0.436 
-1.445 0.442 9/2+ [404] 0.107 


-1.40 0.32 3/2" [541] 0.39 
-1.17 0.33 3/2- [541] 0.33 




^"*Pm 63 
i^^Pm 64 


-0.950 0.422 5/2" [532] 0.789 
-0.783 0.412 5/2" [532] 0.800 


-1.00 0.35 5/2" [532] 0.72 
-0.81 0.35 5/2" [532] 0.74 





Table V. The same as Table H] but for Cs (Z=55) and La {Z—57) isotopes. 





RMF+BCS/TMA 


RHB 


Exp. 


N 


Sp P2 P orbital u"^ 


Sp (32 P orbital ti'' 


Sp 


'''Cs 56 
"^Cs 57 
i^^Cs 58 


-1.913 0.206 1/2+ [420] 0.949 
-1.213 0.171 5/2+ [413] 0.906 
-0.697 0.222 1/2+ [420] 0.94 


-1.97 0.20 l/2+[420] 0.74 
-1.46 0.20 l/2+[420] 0.74 
-0.94 0.21 1/2+ [420] 0.73 


-0.807(7)^^' 
-0.9593(37)***' 


^^'^La 58 
"^La 59 
"^La 59 


-2.136 0.281 1/2" [550] 0.994 
-0.992 0.357 3/2" [541] 0.8 
-0.320 0.343 1/2+ [420] 0.58 


-1.97 0.26 l/2+[420] 0.20 
-1.09 0.30 3/2- [541] 0.73 


-0.806(5)3°' 



energies. Also, it is suggested in Ref.^^^ that such an additional interaction could be 
represented by a surface delta-function interaction. According to our calculations, 
the remaining possible proton emitters in the La isotopic chain are ^^^La and ^^^La, 
as listed in Table V. 

3.6. Quadrupole deformation 

The calculated mass quadrupole deformation parameter, /32m , for odd-Z and 
even-Z nuclei with 54 < Z < 73 at and beyond the proton drip line are plotted in 
Fig. 7 as functions of the neutron number N. While prolate deformations (0.15 < 
/?2m < 0.25) are calculated for Cs isotopes, the proton-rich isotopes of La, Pr, Pm, 
Eu and Tb are strongly deformed (0.35 < (^2 < 0.45). By increasing the number 
of neutrons. Ho isotopes are caused to display a transition from prolate to oblate 
shapes, while most Tm nuclei have oblate shapes. Lu and Ta isotopes exhibit a 
transition from oblate to prolate shapes. The absolute value of /32m decreases as the 
neutron number approaches the conventional magic number, N = 82. 

For stable even-even nuclei, the information concerning the deformation param- 
eter can be derived from measurements of the B{E2) | values. The B{E2) | values 
are basic experimental quantities that do not depend on the nuclear model. Assum- 
ing a uniform charge distribution out to a distance R{9, <j)) and zero charge beyond, 
/?2p is related to B{E2) j by 

/?2p = {4n/3ZRl)[BiE2) T /e^]'/^ (3-4) 

where Rq has been taken to be 1.2A^^^ fm and B{E2) | is in units of e'^b'^. Unfortu- 
nately, we do not have much knowledge about the B{E2) | values for proton drip-line 
nuclei. The two available experimental results^^^ in the region are (32p = 0.221(7) 
for ^^^Xe and /32p = 0.385(48) for ^■^^Ce, which are in reasonable agreement with our 
calculations: /32p = 0.199 for ^^^Xe and /32p = 0.344 for ^"^^Ce (also see Fig. 7). 
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Fig. 7. Ground state quadrupole deformations for proton drip-line nuclei with 54 < Z < 73 as 
functions of the neutron number, A'', at and beyond the proton drip line. 



Another source of knowledge for the deformation of proton drip-hne nuclei can 
be obtained from analysis of the properties of proton emitters, as demonstrated in 
Refs.^^"^^ In those works, to simplify the calculations the assumption that the parent 
and daughter nuclei have the same deformation is used. Our calculations show that 
for Ho, Tm, Lu and Ta isotopes, the parent and daughter nuclei have almost the 
same deformation. However, for the other six isotopes, an absolute deviation of 0.02- 
0.09 is found to exist. To what extent such a deviation can influence the conclusions 
drawn from those calculations^^"^^ needs to be studied in more detail. 



§4. Conclusion 



In the present work, the deformed RMF+BCS model has been used to study the 
ground state properties of the proton drip-line nuclei with 55 < Z < 73, including 
the location of the proton drip line, the ground state quadrupole deformation, the 
one-proton separation energy, the deformed single-particle orbital occupied by the 
odd valence proton, and the corresponding spectroscopic factor. The RMF+BCS 
model reproduces the available experimental data reasonably well, except in the 
case of odd-odd nuclei. The results also agree well with those of the RHB method. 
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The systematic discrepancies for the one-proton separation energies of the odd-odd 
proton-rich nuclei may be attributed to the lack of a residual proton-neutron pairing 
in both mean field models. It will be very interesting to include a residual proton- 
neutron pairing into the relativistic mean field model through a generalized BCS 
method in a future work. 

In conclusion, we have found that the RMF-I-BCS model can describe the ground 
state properties of proton emitters as well as the more complicated RHB model. It 
is also found that the state-dependent BCS method with a zero-range (5-force is valid 
not only in the stable region but also in the proton drip line. To summarize, the drip- 
line nuclei predicted by the RMF-^BCS/TMA calculations are ^^^Ta, ^^^Lu, ^^^Tm, 
i44jjo, 139Tb, i^^Eu, i29pni, i25pj.^ ii8La and ^^^Cs. Based on our calculations, 
possible proton emitters that could be detected in future experiments are ^^^Lu, 
^^^Lu, ^^^Lu and ^^^Lu for the Lu isotopic chain, ^'^^Ta for the Ta isotopic chain, 
^'^°Eu for the Eu isotopic chain, ^^^Tb and ^^^Tb for the Tb isotopic chain, i^^Pr and 
I20pj. ^Yie Pr isotopic chain, ^^"^Pm and ^^^Pm for the Pm isotopic chain, -'■^■'■Cs 
for the Cs isotopic chain, and ^^^La and ^^^La for the La isotopic chain. 
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